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Application of an Elongated Kelvin 
Model to Space Shuttle Foams 

Spray-on foam insulation is applied to the exterior of the Space Shuttle’s External Tank to limit 
propellant boil-off and to prevent ice formation. The Space Shuttle foams are rigid closed-cell 
polyurethane foams. The two foams used most extensively on the Space Shuttle External Tank are 
BX-265 and NCFI24-124. Since the catastrophic loss of the Space Shuttle Columbia, numerous 
studies have been conducted to mitigate the likelihood and the severity of foam shedding during the 
Shuttle’s ascent to space. Due to the foaming and rising process, the foam microstructures are 
elongated in the rise direction. As a result, these two foams exhibit a non-isotropic mechanical 
behavior. In this paper, a detailed microstructural characterization of the two foams is presented. 
The key features of the foam cells are summarized and the average cell dimensions in the two 
foams are compared. Experimental studies to measure the room temperature mechanical response 
of the two foams in the two principal material directions (parallel to the rise and perpendicular to the 
rise) are also reported. The measured elastic modulus, proportional limit stress, ultimate tensile 
stress and the Poisson’s ratios for the two foams are compared. The generalized elongated Kelvin 
foam model previously developed by the authors is reviewed and the equations which result from 
this model are presented. The resulting equations show that the ratio of the elastic modulus in the 
rise direction to that in the perpendicular-to-rise direction as well as the ratio of the strengths in the 
two material directions is only a function of the microstructural dimensions. Using the measured 
microstructural dimensions and the measured stiffness ratio, the foam tensile strength ratio and 
Poisson’s ratios are predicted for both foams. The predicted tensile strength ratio is in close 
agreement with the measured strength ratios for both BX-265 and NCFI24-124. The comparison 
between the predicted Poisson’s ratios and the measured values is not as favorable. 
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Acreage foam is NCFI24-124 
(Machine sprayed) 


Close-outs 
BX-265 or 


are hand sprayed 
PDL-1034 
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Typical Flight Loads on ET Foam Applications 
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Foam Microstructure 



• 97% air; y = 0.03 

• polymeric cell walls 

• due to its micro structure, material is anisotropic (possess 
different material properties in different directions) 



Parallel-to-rise 

direction 


> 

Perpendicular-to-rise 

direction 


Figure 3.5 - Cell Geometry, NCFI24-124 
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The foam microstructure can be approximated by an elongated 
tetrakaidecahedron (a fourteen-sided polyhedron) 




NCFI24-124 Average number of 

faces per cell: 13.7* 


* Wright L. S. and Lerch B. A., 2005. Characterization of space shuttle insulative materials, NASA/TM-2005-2 13596. 
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Some Significant Previous Studies 

The tetrakaidecahedron foam model is commonly referred to as the Kelvin foam model after: 

Thomson, W. (Lord Kelvin), 1887. On the division of space with minimum partitional 
area. Phil. Mag. 24, 503-514. 



William Thomson (Lord Kelvin) determined that the tetrakaidecahedron (with slightly 
curved faces) was nature ’s preferred shape for soap bubbles and other foams since it is 
the shape that minimizes the surface area per unit volume and packs to fill space. 


Equi-axial Unit 
Cell Models 


Zhu H.X., Knott J.F. and Mills N.J., 1997. Analysis of the elastic properties of open-cell foams with 

tetrakaidecahedral cells. Journal of the Mechanics and Physics of Solids 45 (3), 319-343. 
Warren, W.E. and Kraynik, A.M., 1997. Linear elastic behavior of a low-density Kelvin foam with 
open cells. Journal of Applied Mechanics 64 , 787-794. 


Elongated Unit 
Cell Models 


r 

< 




Dement’ev, A.G. and Tarakanov, O.G., 1970. Model analysis of the cellular structure of plastic foams of the 

polyurethane type. Mekhanika Polimerov 5, 859-865. translation Polymer Mechanics 6, 744-749. 

Gong, L., Kyriakides, S. and Triantafyllidis, N., 2005. On the stability of Kelvin cell foams under compressive 
loads. Journal of the Mechanics and Physics of Solids 53, 771-794. 

Ridha, M., Shim,V.P.W. and Yang, L.M., 2006. An elongated tetrakaidecahedral cell model for fracture in rigid 
polyurethane foam. Key Engineering Materials 306 - 308 , 43-48. 
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The size and shape of an elongated tetrakaidecahedron are uniquely defined 
by specifying three independent dimensions. 



H = 4L sin 6 
D = 2L cos 6 + -Jib 


Aspect Ratio of Unit Cell 

R _H _ 4L sin <9 
D 2 L cos 6 + fib 
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The majority of the solid mass resides in the cell edges (where the faces 
come together). 




Face thickness in the middle of the faces 
~ 0.1 jam to 1.0 jam 


Simplifying Assumptions: 

1) The structural rigidity of the cell faces is assumed to contribute little to the foam 
mechanical behavior. 

2) The mass of the cell faces are only a small fraction of the total solid mass. 
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The mechanical behavior can be accurately modeled by considering the deformation of 
the cell edges only. The cell edges are assumed to act like struts possessing axial and 
flexural rigidity. 


Z (rise direction) 






Relative Density 

P Km 8A(2 L + b) 242 L + b) 

Ps Kota! HD £ sin 0 {lL COS 0 + 'flhf 
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Cell edge cross-sections are approximated as three-cusp hypocycloids (Plateau Borders) 



Cross-section properties 

lx ■> ly bending moments of inertia 
A Cross-sectional area 
S X ,S Section moduli 

are a function of the cross-section radius only. 



A = {^[3 -7r/2jr 2 
I x -I y - (20a/3 - 1 \7rjr 4 / 24 
S x = (60 - 1 1 V3;r) r 3 / 24 
S y =(20V3-ll^-)r 3 /12 


Thus, four microstructural dimensions are required: 

3 to define the unit cell + 1 to specify the edge cross-section dimension 
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Careful consideration of the unit cell deformation leads to convenient algebraic 
equations for the Young ’s modulus, Poisson *s ratio and strengths in the principal 
material directions. 



Poisson’s ratios: 


b[Ab 2 -\2l) 

~~ Vyx ~ 12/(21 cos 2 ~9 + b]+ ffp sin 2 <9 + b 3 ) 


Young’s moduli: 


C=Ey=- 


12 El 


/sin# 


2/ 3 sin 2 0 + b 3 H ( 2 / cos 

A V 


E, = 


24// sin 0 


2 ^ 12/ sin 2 # 

cos # + — 

AL 2 


[ 12 . 


2/cos# + 6 


Strengths: 


r 

based on peak stress in any edge 
reaching solid material strength G s 

v J 


yy 


L cos 6 sin 6 


- + - 


Z 2 sin 2 6 

255 


2Zcos# + 


V2 b] 


U xz ~ ^yz ~ ' 


(aP - \2I%L cos # + f2b\ cos 0 


12/(2/ cos 2 0 + b)+ A(2L 3 sin 2 0 + b 3 


a 


cr v 


<j 


yy 


L sin 0 Lb sin 6 
+ 


flA 2V25 


2L cos 0 + 


V2 b] 


ffP (AL 2 - 1 2/)cos # sin 2 6 


[l2/Lsin 2 0 + AL 3 cos 2 #] 

dLcosO + b 

® zz 

sin 6 L cos # 

\flL cos # + b^ 




| 

2 A 45* 




e+b) 


Failure of the 
L-length edges 


Failure of the 
b-length edges 
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The ratio of the rise direction modulus to the normal-to-rise direction modulus 
and the ratio of the strengths are indicative of the amount of elongation. 



Stiffness Ratio 



E z _R 2 

2 sin 2 0 + (b/ L) 3 cos 2 

L al 2X 

0 + b/L ) 

Ey 4 

2 n 12/ . 2 n 

cos 6 + — z-sin 6 

l aT J 




Strength Ratio 


Doesn ’t involve the properties of the solid material. 
Only a function of the unit cell dimensions and 
the edge cross-section dimension. 


° s zz ° s zz D cos 0 + AL sin 6 

<J ~ _ — K 7 

<j s xx a S yy 2 S x sin O + AL cos 6 
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Microstructural Characterization 






Cell Height H 

Cell Width D 

Edge Cross- 
section Radius r 

Aspect 
Ratio 
R = H/D 

Avg., 

|Lxm 

n 

Avg., 

|um 

n 

Avg., 

|Lim 

n 

BX-265 

193 

100 

136 

100 

NA* 


1.42 

NCFI24-124 

248 

100 

142 

100 

26.0 

27 

1.75 


* Data not available 
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Mechanical Testing 


BX-265 





Average Properties 



Initial Modulus 
(MPa) 

Proportional Limit 
(kPa) 

Ult. Tensile Strength 
(kPa) 

Ratios 

Rise 

Normal 

-to-rise 

Rise 

Normal- 

to-rise 

Rise 

Normal- 

to-rise 

Stiffness 

Proportional 

Limit 

Ultimate 

Tensile 

Strength 

BX-265 

13.53 

7.03 

315.4 

173.5 

555.8 

319.4 

1.92 

1.82 

1.74 

NCFI24- 

124 

20.80 

7.07 

248.1 

110.3 

375.4 

188.2 

2.94 

2.25 

1.99 
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Mechanical Testing 



Poisson ’s Ratios 


Average Density 



Average 

Density 

(g/cc) 

Relative 

Density* 

BX-265 

0.0369 

0.031 

NCFI24-124 

0.0373 

0.031 


* Calculated assuming density of polymer 
is 1.2 g/cc. 




Avg. Measured 

Standard 



Poisson’s Ratio 

Deviation 


V xy 

0.355 

0.06 

BX-265 

V xz 

0.273 

0.0007 


V zx 

0.536 

0.29 



0.675 

0.17 


V xy 

0.382 

0.14 

NCFI24-124 

V xz 

0.183 

0.04 


V zx 

0.641 

0.10 
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The shape is defined by the aspect ratio 


and the shape parameter Q = 

L cos# 



Ra > Rb 





Qa > Qb 
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in terms of L, b, 6, I/A 


r 4 L sin 6 ^ 

2sin 2 9 + (b/Lf+ UI 2 (2 cos 2 9 + b/ l) 
_ A.L _ 

K 2Lcos9 + /2b 2 

4 

2/1 12/ -2/1 
cos 9 + — z- sin 9 

AL 2 




R = 


4 L sin 6 
2 L cos 6 + Vl b 


Q 


L cos# 


in terms of R, Q, y 



r e = — 

E 4 


2 Q 2 R 2 + 64 ^ 


c | 8i?C 2 g 3 (32 + 4gVl6 + g 2 i? 2 / ^ 


Vi 6 +e 2 i ? 2 J ( 4 e+ 2 Vi 6 +e 2 J R 2 )(i 6 +e 2 J R 2 ) 


16 Cj + 


8R 3 C 2 Q 5 


4Q + 2 A /l6 + (? 2 ^ 2 ) (l 6 + g 2 ^ 2 ) 
Q = 2 + /2Q 


r 


Q = 2 + f>Q 



^ 4 L sin $ ' 

2L cos 0 + 4lb 


f Sf 

AL 


r 2Sf 

AL 


\ 

cos 9+ sin 9 

l 

\ 

sin 0 + cos 6 

J 


r = 


2A(2L + b) 


L sin 9 ( 2 L cos 9 + J2bj 


OR A 

16V2C 3 g 1 - 5 i? 0 V 0 ' 5 


(4Q + 2 f 6 + Q 2 K 2 f 5 (l6 + Q 2 R 2 Jf 

4 A /Ci + 

4V2C 3 g 2 - 5 i? L V 0 ' 5 


(4e+2Vi6+e 2 i? 2 ) )j (i6+e 2 i? 2 ) 0 ' 5 _ 
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Previous research studies using the elongated Kelvin model 
restriction on the unit cell shape, 


b! L = V2 cos 6 



q=4i 


which is equivalent to assuming that 


D = 24lb 


This reduces the number of unit cell dimensions 
required to exercise the theory and apply the equations 
by one. The unit cell shape is now described by 
specifying only two dimensions. 


set an artificial 




However, this restriction on the unit cell geometry seems arbitrary and it reduces the 
generality of the model and limits its applicability to a narrower range of foams. 
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Plateau Border Cross-section 




8 


6 

LU 

0 5 

-i— > 

CO 

“ 4 

CO H 

(Ji 

CD 

1 3 

CO 


0 


□ Gupta et al. (1986), Polyisocyanurate 
OHilyard (1982), Flexible Polyurethane 
A Huber & Gibson (1988), Flexible Polyurethane 
O Huber & Gibson (1988), Rigid Polyurethane 

♦ Mehta & Colombo (1976), Polystyrene 
XGong et al. (2005a), Polyester-Urethane 
▲ BX-265 Rigid Polyurethane 

♦ NCFI 24-124 Rigid Polyurethane 



1.0 


1.2 


1.4 1.6 

Aspect Ratio R 


1.8 


2.0 
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D(jum) 

R 

r 

Re 

BX-265 

136 

1.42 

0.031 

1.9 

NCFI24-124 

142 

1.75 

0.031 

2.9 



l 




b = 


D 


4i + - 

Q 


L = 


_ b^l 16 + (2 + ^qJ R : 

4 Q 


r = 


y RD : 


C x (l 6L + 8Z?) 



BX-265 

NCFI24-124 


Predicted 

Measured 

Predicted 

Measured 

Q 

1 .0747 

— 

0.7765 

~ 

K 

1.75 

1 .82 PL 
1.74 ult. 

2.34 

2.25 PL 
1.99 ult. 

b {/um) 

41.5 

35.0* 

35.6 

NA 

L (/an) 

61.8 

63.0* 

77.2 

NA 

r (/an) 

22.8 

18.0* 

25.2 

26.0 

»xy 

.177 

.35510.06 

.060 

.38210.14 


.422 

.27310.0007 

.373 

.18310.04 


.811 

.53610.29 

.67510.17 

1.097 

.64110.10 


* Measurements from another block of BX-265. 
NA Data not available. 
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Finite Element Analysis of an Elongated Kelvin Model 

Unit cell dimensions 




Closed Cell 



R = 1.72 
D = 144 jum 
b = 24 jLim 
L = 90.6 jum 
r = 20.2 (jm 

Solid Properties 

E = 2GPa 
v = 0.3 



Open Cell 



Closed cell 

Open cell 

y 

1 .62% 

1 .24% 

V X y 

0.621 

0.024 

N 

X 

> 

0.120 

0.230 

v 

v zx 

0.256 

1.727 
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Summary and Conclusions 



• Model was successful in predicting the measured strength ratios and edge cross-section 

radii using measured average cell height, cell width, relative density and stiffness 
ratios as input. 

• Prediction of Poisson’s ratios was not as successful. 

• Shuttle foams have a microstructure such that Q * V 2 ? so a micromechanics model 

derived from a general elongated Kelvin unit cell is needed. 

• Future work should include performing finite element analysis of BX-265 and 

NCFI24-124 unit cells with faces included and predict Poisson’s ratios. 
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Backup Slides 
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Manual Spray Close-outs of BX-265 



Intertank Flanges 



Protrusion Air Load (PAL) Ramp 


Redesign eliminated the LOx 
and LH 2 PAL Ramps 
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Current design 
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Bipod Fitting and 
Ramp Closeout 

BX-265 
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Ice Frost Ramps 

PDL-1034 

Cover the pressurization line and cable tray 
support fittings to prevent ice formation. 
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Structural Analysis of the Foam Applications under launch loads and environments 
is performed by Finite Element Analysis 





Ej 

PSI 

m _ 

r 111300 

16.142 

w 

l 101600 

14.735 


- 31040 

13.320 


- 02030 

11.306 


- 72350 

10.433 


- 62610 

3.000 


- 52070 

7.660 


- 43120 

6.254 


- 33300 

4.041 


- 23640 

3.423 


- 13030 

2.015 


4152 

0.602 


L -5531 

-0.011 



Finite element analysis treats the foam material as a homogeneous, temperature dependant 
and orthotropic material 

• elastic constants and strengths obtained during material characterization testing 

• the effect of a varying micro structure on the foam structural integrity is neglected 


Analyses are not used for flight qualification, but as a tool to study possible foam 
shedding mechanisms and to guide proposed design changes. 
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